FINE REGULARITY FOR ELLIPTIC SYSTEMS 
WITH DISCONTINUOUS INGREDIENTS 

DIAN PALAGACHEV AND LUBOMIRA SOFTOVA 

Abstract. We propose results on interior Morrey, BMO and Holder regularity for the 
strong solutions to linear elliptic systems of order 2b with discontinuous coefficients and 
right-hand sides belonging to the Morrey space L P ' X . 



1. Main Results 

It is well known that, when dealing with elliptic systems, in contrast to the case of 
a single second-order elliptic equation, the solely essential boundedness of the principal 
coefficients is not sufficient to ensure Holder continuity even of the solution (see [HI Chap- 
ter 1]). On the other hand, precise estimates on Holder's seminorms of the solution and 
its lower order derivatives is a matter of great concern in the study of nonlinear elliptic 
systems. In fact, these bounds imply good mapping properties of certain Caratheodory 
operators ensuring this way the possibility to apply the powerful tools of the nonlinear 
analysis and differential calculus. It turns out that "suitable continuity" of the principal 
coefficients of the system under consideration is sufficient to guarantee good regularity 
(e.g. Sobolev) of the solutions (see [5, 12J). We deal here with discontinuous coefficients 
systems for which the discontinuity is expressed in terms of appurtenance to the class 
of functions with vanishing mean oscillation. Although such systems have been already 
studied in Sobolev spaces W 2b ' p (cf. [3J), our functional framework is that of the Sobolev- 
Morrey classes W 2b,p,x . These possess better embedding properties into Holder spaces than 
W 2b,p and as outgrowth of suitable Caccioppoli-type estimates, we give precise character- 
ization of the Morrey, BMO or Holder regularity of the solution and its derivatives up to 
order 2b — 1 . 

Let Q be a domain in M. n , n > 2, and consider the linear system 
(1.1) £(x, D)u := J2 Aa(x)D a u(x) = f (x) 

\a\=2b 

for the unknown vector- valued function u: Q — ► IR m given by the transpose u(x) = 
(ui(x), . . . , u m (x)) , m > 1, f(x) = (fx(x), . . . , f m (x)) T , where A a (x) is the m x m- 
matrix {a^ fc (x)}™ fc _ 1 and a?^\ Vt — > K. are measurable functions. Hereafter, b > 1 is 
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a fixed integer, a = (ax, . . . ,a n ) is a multiindex of length |a| = a± + ■ ■ ■ + a n and 
D a := D" 1 . . . D% n with Di := d/dxi. This way, the matrix differential operator £(x, D) 
has entries 

P k (x,D):= a a( x ) Da 

\a\=2b 

and for fixed j and k the polynomial 

|a|=26 

is homogeneous of degree 2b. 

We suppose (I1.1J) to be an elliptic system, that is, the characteristic determinant of 
£(x, £) is non-vanishing for a. a. x E Q and all £ 7^ 0. In view of the homogeneity of P fc 's, 
this rewrites as (see |H1E|) 

(1.2) 3 5 > 0: det { ^ A Q (x)c} > S\^\ 2brn a.a. V£ G M n . 

|a|=26 

Our goal is to obtain interior Holder regularity of the solutions to (|1.1|) as a byproduct 
of a'priori estimates in Sobolev and Sobolev-Morrey spaces. Let us recall the definitions 
of these functional classes. 

Definition 1.1. The Sobolev space W 2b,p (fl), p G (l,+oo), is the collection of L P (Q) 
functions u: Q — > R all of which distribution derivatives D a u with \a\ < 2b, belong to 
L P (Q). The norm in W 2Wp (Q) is 

2b . v l/p 

\\u\\w^p(n) '■= ^ \\D a u\\p ; n, \\ ■ \\ p; n := ( / | • \ p dx J 

s=0 \a\=s ^ Q ' 

For the sake of brevity, the cross-product of m copies of L P (Q) is denoted by the same 
symbol. Thus, if u — (u 1 ,...,u m ) is a vector-valued function, u G L P (Q) means that 
u k G L p (Vt) for all k = 1, . . . ,m, and \\u\\ p . n := YJk=i IKIUn- 

Definition 1.2. Let p G (l,+oo) and A G (0, n). The function u G L P (Q) belongs to the 
Morrey space L P ' X (Q) if 

[If y/ p 

SUp — r- / |M(x)| P (ix < OO 

\r>o r A J Brnn J 



\ u \\p,\;n 



,r>o r JB r nn 

where B r ranges in the set of balls with radius r in M. n . The Sobolev-Morrey space 
W 2b,p,x (Q) consists of all functions u G W 2b,p (Q) with generalized derivatives D a u, \a\ < 
2b, belonging to L P,X (Q). The norm in W 2b ' p ' x (Q) is given by 

2b 

\\u\\ W 2b, P ,\(ji) := ^ ||-D a M||p,A;f2- 
s=0 \a\=s 
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We refer the reader to El Ej for various properties of the Morrey and Sobolev-Morrey 
spaces. 

Definition 1.3. For a locally integrable function /: W n —>■ R define 

r]f{R) : = sup — !— / \f(y) - f Br \dy for every R > 0, 

r<R \-t3 r \ J B r 

where B T ranges over the balls in W 1 and fs r = j-J-y j B f{y)dy. Then: 

• / G BMO (bounded mean oscillation, see John-Nirenberg [7j) if ||/||* := sup R r]f(R) < 
+oo. Il/H* is a norm in BMO modulo constant functions under which BMO is a 
Banach space. 

• / G VMO (vanishing mean oscillation, see Sarason ^T]) if / G BMO and 
liniRjo ^/(-R) — 0. The quantity r]f(R) is referred to as KMO-modulus of /. 

The spaces BMO(Q) and VMO(Q), and || • are defined in a similar manner taking 
B r fl Q instead of B r above. 

As already mentioned, the desired Holder bounds will be derived on the base of the 
following a'priori estimate in Sobolev-Morrey classes. 

Theorem 1.4. Suppose a{ k G VMO(tt) n L°°(tt), f G 1 < p < oo, < 

A < n, and let u G W^ P,X (Q) be a strong solution of T/ien, /or any Q' <<= Q" <e 

i/iere zs a constant C depending on n, m, b, p, A, 6, Wa^W^n, the V MO -moduli rj aJ k of 
the coefficients (cf. 011]) and dist (fi',<9f2"), such that 

(1-3) Hull^a&.p^Q/) < C (||f ||p,A ; n" + || u ||p,A ; n") • 

It turns out, moreover, that the operator £ improves the integrability of solutions to 
(jl.lj) . In fact, by means of standard homotopy arguments and making use of formula 
(EH) (cf. 0, PHI Section 3]), it is easy to get 

Corollary 1.5. Under the hypotheses of Theorem \l-4i suppose u G W^! q (Q) with q G 
(l,p). ThenueW^ p >\n). 

A combination of (|1.3|) with the embedding properties of Sobolev-Morrey spaces leads 
to a precise characterization of the Morrey, BMO and Holder regularity of the solutions 

to (JEE). 

Corollary 1.6. Under the hypotheses of Theorem \l.J\ define so as the least non-negative 
integer such that 2b ™ so > 1 and fix an s G {s , . . . , 26 — 1}. Then there is a constant C 
such that: 

a) ifp G (l, f^) then D s u G L^-^+^O') and 

||-D s u||p,(2&-s)p+A ; n' < C (||f||p,A ; n" + ||u||p,A;n") ; 
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b) i/p = ^ ffren ZPu e BMO{VL') and 

\\D S \l\ [ BMO;VL' < C (||f||p,A;n" + || u ||p,A;^") i 

c) »/pe(g, ^^i)* i/ien J ° Su G C° ,<Ts ( fi/ ) «»'^ cr s = 26 - s - ^ and 




26 - s " ' 26 - s 26 - s - 1 



The picture illustrates geometrically the results of Corollary 1 1 . 61 with the couple (p, A) 
lying in the semistrip |(p, A) : p > 1, < A < n} and s G {so, . . . , 26 — 1}. The points B s 
on the p-axis are B s = (^r, 0) , 5 = (1, 0), and A s = (1, n — 26 + s) is the intersection 
of the line through (0, n) and B s with the vertical line {p = 1}. 

When (p, A) belongs to the open right triangle BB S A S then we have case a), i.e. 
D s u E L p >( 2b - S 1 p+X (tt'). If (p, A) lies on the open segment A S B S we have D s u E BMO(Q') 



This rewrites asp £ (n — A, oo) when s = 26 — 1. 
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(case b)). In particular, (p, A) G ABB S0 A S0 yields L> So u G LP'^- 30 ^^'), while u G 

if Sq > i Furt her on, (p, A) G A so £ So gives L> s °u G BMO{VL'). 

Let s G {s , • • • , 2b — 2} and suppose (p, A) lies in the interior of the quadrilateral 
Q s := B S B S+1 A S+1 A S . Then D s u G C°' a °(tt') (case c)) whereas D s+1 u G L^ 26 -*- 1 )^^'). 
Moreover, the exponent a s is the length of the segment CA S where C = C(p,X) is the 
intersection of the line {p = 1} with the line passing through the points (p, A) and (0, n). 
In particular, (p, A) G A s+1j B s+1 implies £> s+1 u G BMO(Vt'). 

Similarly, set Qib-\ for the shadowed unbounded region on the picture. Then (p', A') G 
Q 26 _i gives -D 26 " 1 !! G L7°- CT ^-i(fi') with a 2b -i equals to the length of C'A 2b -u C = 
C'(p', A'), while D 2b - l u G BMO{Q!) if (p', A') G A 26 _ lJ B 26 _ 1 . 



2. Newtonian- type Potentials 

Fix the coefficients of (jl.ljl at a point s G and consider the constant coefficients 
operator £(xo,D) := J2\ a \=2b-^a(xo)D a . The 26m-order differential operator 

L(x ,D) := det£(a;o,£>) = det | ^ a^(x )D a ^ 

\a\=2b 

is elliptic in view of (jl.2J) . and therefore there exists its fundamental solution T(xq; x — y) 
(see jSJE]). If the space dimension n is an odd number, then 

f (x ; x-y) = \x- y\ 2bm - n P (x ; 

V \x-y\ 

where P(xq]^) is a real analytic function of £ G § n_1 := {£ G K": |£| = 1}. In case of 
even dimension n the procedure is standard and is based on introduction of a fictitious 
new variable x n+ i and extension of all functions as constants with respect to it, see jH] for 
details. Set {Lj k (x , £)}™ fc=1 for the cofactor matrix of {P k (x Q , £)}Jj. =1 and note that for 
all fixed j,k = l,...,m, Lj^(xq, D) is either a differential operator of order 2b[m — 1), or 
the operator of multiplication by 0. Making use of the identities 

rn 

fe=l 

with Kronecker's 5ij, it is not hard to check (cf. [5J that the fundamental matrix 
T(xo;x) = {P fc (xo; x)}™ fc=1 of £(xo,D) has entries P fc (xo;x) = Lkj(xo, D)T(xq; x). 

Let r > be so small that B r = {x G R n : \x - x \ < r} m tt, and let v G C£°(B r ). 
Then, employing 



£(x , D)v(x) = (£(x , D) - £(x, D))v(x) + £(x, D)v(x) 
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and using standard approach (cf. 0E|), we obtain a representation of v in terms of the 
Newtonian-type potentials 



via; 



/ T(x ;x-y)£v(y)dy+ f T(x ; x — y) (£(x , D) — £(y, D))v(y)dy. 

J Bf J Br 



Taking the 26-order derivatives and then unfreezing the coefficients by putting x = x, 
we get 

(2.1) D a v(x)=p.v. D a T{x;x -y)£v{y)dy 

J Br 

V v ' 

=: #a(£v) 

+ P v - [ D a T(x;x-y)(A a ,(x)-A a/ (y))D a 'v(y)dy 

=■ e a [A Q ,,D a 'v] 
+ / D pa r(x;y)v s da y £v(x) Va:|a|=26 

where the derivatives D a T(-; •) are taken with respect to the second variable, the multi- 
indices /3 s are such that (3 s := (aj, . . . , a s _i, ot s — 1, a s+ i, . . . , a n ), \P S \ = 26 — 1 and u s is 
the s-th component of the outer normal to S n_1 . 

Noting that each entry of the matrix D a T(x; y), \a\ = 2b, is a C alder on- Zygmund kernel 
(cf. OIU), we have 

Lemma 2.1. Let \ot\ = \a'\ = 2b and A a G L°°(Q). For each p G (l,oo) and each 
A G (0, n) t/iere a constant C depending on n, m, b, 5, HA^H^n, p and A such that 

(2.2) ||^af \\p,X;n < C||f||p,A;^' ||^a[-A- Q ', f] ||p,A;fi < C||A Q / ||* ; n||f ||p,A;fi 

/or a// f G L p ' A (fi). Moreover, Zei A a G MO(O) nL°°(fi) iw't/i VMO-modulus r} Aa . Then 
for each e > i/iere exists r = ro(e,r)A a ) such that if r < r we have 

(2-3) \\€ a [A a ,,i]\\ P ,X;B r <Ce\\i\\ pXBr 

for all B r m Q. and all f G L p ' x (B r ). 

Lemma l2~T1 is proved in ^UJ Theorem 2.1, Corollary 2.7] in the general case of Calderon- 
Zygmund's kernels of mixed homogeneity. 

3. Proof of Theorem 11.41 

Fix an arbitrary xq G suppu and let B r := {x G M™ : \x — Xq\ < r}. The operators 
and (£ a [A a ', •] are bounded from LP into itself (cf. [3]) and therefore the representation 
(|2.1|) still holds true (a.e. in Q) for any function v G W ' p (fi) := closure W 2b,p^C^°(Q) 
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and moreover for v G WQ b,p (fl) H W^ff ' (ft) as well. Let supp v C B r . In view of (j2.1j) . 
()2.2|) and ()2.3|) . for each £ > there exists r (e, r/A a ) such that for r < r one has 

\\D 2 \\\ p>x , Br < C(\\2v\\ p , x . Br + e\\D 2b v\\ p , x . Br ) 

whence, choosing e small enough we obtain 

(3-1) \\D 2b v\\ M , < C\\2v\\ PjX]Br . 

Let G (0,1), 9' = 0(3 - 0)/2 > and define a cut-off function (p(x) G C^{B r ) such 
that </?(x) = 1 for x G -Bgr whereas v ? ( a; ) = for x G" Bqi t . It is clear that \D s (p\ < 
C(s)[0(l - 9)r]- s for any 1 < s < 2b because of 9' - 9 = 9(1 - 0)/2. Applying Jjjj to 
v(z) := y?(x)u(x) G Wq 26 '^) n W 26 'f- A (S r ), we get 

||-D 26u ||p,A;B flr < \\D 2b v\\ P: x- Bg , r < C\\£,v\\ P: \.B g , r 

(26—1 I. n 26— s II II II 

|| f || \- 11-^ U llp,A;B flV l|u||p,A;B 9V 

||1 \\p,\;B glr f ^ ^ _ f _ 

Hence, the choice of 0' and 0(1 - 0) < 20'(1 - 0') imply 

[0(1 - e)r] 2b \\D 2b u\\ p>X;Ber < c([0'(l - 0')r] 2 "||f || p , A;BeV 

26-1 

+ - / )r] s ||L> s u|| PjA . Be , r + ||u|| PiA;BeV ) 



s=l 



Setting 6 S for the weighted Morrey seminorms sup o<e<1 [0(l — 9)r] s \\D s u\\ PtX - Bf)r , s G 
{0, . . . , 26}, the last inequality rewrites as 



26-1 



(3.2) e 2b < c ^r 26 ||f || p , A;Br + ®s + e J . 

To manage the seminorms G s with 1 < s < 2b — 1 , we use the next interpolation inequality 
which follows from [T2*l (5.6)] as in [TTH Proposition 3.2]: 

Proposition 3.1. There exists a constant C = C(n,m,b,p, A, s) independent of r and 
such that 

e s < e& 2b + -^Lye VeG(0,2). 

Therefore, interpolating the intermediate seminorms in ()3.2|) and fixing = 1/2, we 
obtain the following Caccioppoli-type inequality 

(3-3) p 26 u|| PjA;Br/2 < C(||f|| PiA;Br . +r- 26 ||u|| PiA;Br ). 

The estimate ()1.3|) follows from ()3.3|) by a finite covering of Q' with balls B r /2, r < 
dist(fi',0ft"). 
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4. Proof of Corollary II .HI 

To begin with, let s = 2b - 1 and fij. = 5 r n ft' with 2r < dist (fi', dQ"). Direct 
calculations based on the classical Poincare inequality lead to 

(4.1) JL / |D^y a ;)-(D^ 1 u)n;|^<C(n,i» I m)||D a6 u||J iAsgl 

with suitable f2' such that Q' d fi' <s fi" and (f)n^ standing for the integral average of 
f over Q' r . Taking the supremum with respect to r we get that D 2h ~ l u belongs to the 
Campanato space £ P ' P+A (f2'). Now, employing the embedding properties of Campanato 
spaces into Morrey, BMO and Holder ones (cf. QjJ Theorem 2.1]) and (jl.l-jj) . we obtain as 
follows. If p + X < n then D 2h ~ l u G L p ' p+ \tt')^ p + X = n implies D 2h ~ l u G BMO (ft'), 
while D 26_1 u G C°' CT2b - 1 (fi 7 ) with (T 2 6-i = 1 - (n - A)/p when p + A > n. 

To prove the statement for any s, we run induction for decreasing values of s until sq. 
Thus, suppose D s+l u satisfies the statement of Corollary 11.61 Then 



, — Y~rr f \D s u(x) - (D s u) n ,\ p dx < \\D s+1 u\\ p r s - 

r (26-s)p+A J I V / V ^ l r\ — M ll p,(2b-s-l)p+X;n' 



<^(l|f|liw + IWw) 

whence D s u G £ P '( 26_S ' >P+A (r2 / ) and the conclusions follow as above. 

Similarly, if s > 1 and p G (l, ^E^) we have D^u G £p-(26-so+i)p+A(q/)_ It fol _ 
lows from the definition of Sq that n < (2b — so + l)p + A < n + p, whence D s °~ 1 u G 

c°> 2b - so+1 -^(n'). 
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